Correlation length for amorphous systems 
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Crystals and quasicrystals can be characterized by an order that is a purely geometric property 
of an instantaneous configuration, independent of particle dynamics or interactions. Glasses, on 
the other hand, are ostensibly amorphous arrangements of particles. A natural and long-standing 
question has been whether they too have, albeit in a hidden way, some form of geometric order. Here 
we define a coherence length that applies to systems which are typically characterized as amorphous, 
as well as to those that are conventionally ordered. We argue that the divergence of such a length 
is consistent with current theories of the 'ideal' glass transition. 
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There comes a moment in many introductory talks on 
glasses when two pictures are projected showing config- 
urations of a liquid and a glass, and the speaker defies 
the audience to tell which is which. The purpose of this 
paper is to propose a purely geometric correlation length 
that may indeed do precisely this, characterizing glass 
structure even though we may not know, a priori, what 
the ideal glass state is. The idea is to consider the spa- 
tial frequency with which any given patch of linear size 
£ and volume V ^ appears in an infinite sample. We 
shall argue in what follows that patches smaller than a 
certain size io occur 'often' - with an average repeat dis- 
tance V that scales sub-exponentially with their volume, 
while patches larger than £o do so 'rarely', with V being 
exponential in V, i.e., just as one would expect of a ran- 
dom system. The existence of an 'ideal glass state' then 
involves a divergence of this crossover length £o (with a 
concomitant diverging time scale) below the thermody- 
namic transition temperature T^, implying that in such 
a state motifs of all sizes recur frequently. We note that 
our definition makes no reference to either interactions or 
equilibration, nor is it restricted to systems of a specific 
type of constituent {e.g., particle, polymer, or continu- 
ous density). In the rest of the paper we shall proceed, 
in the context of several examples, in two steps: first, we 
argue for the existence of an ordered io = oo phase where 
all patterns recur frequently, and second, we describe a 
mechanism whereby this length becomes finite. 

That some form of amorphous order exists in the glass 
phase is implicit in several theories of fragile glasses, most 
notably in the 'Random First Order' scenario [11 [J. One 
assumes that there is some 'ideal glass' state to which the 
system would equilibrate, given an infinitely slow anneal- 
ing from high temperature to T < T^, where Tk is the 
Kauzmann temperature. Such an ideal state should be 
characterised by permanent non-periodic spatial density- 
modulations. In this ideal glass, the entropy, or, more 
precisely, the complexity (the logarithm of the number 
of average density profiles) [3^ is subextensive in the size 
of the system. The main point of this paper is that the 



further assumption that in the ideal glass state the com- 
plexity of patches of volume V which may be found in an 
infinite system [4 is subextensive may be used to define a 
purely geometric length scale. To illustrate this point, we 
shall discuss, in some detail, a more stringent variant of 
this constraint, in which the state of a macroscopic por- 
tion of a sample determines, if not completely, at least 
significantly, the adjacent region. 

A thought experiment illustrating this [SJ [H [71 [8] pro- 
ceeds as follows: First, equilibrate a large sample and 
choose a region S of linear size £ and volume V (x De- 
note the (microscopic) configuration of this region by A, 
and look at another region S', contiguous with H, whose 
volume V = aV is proportional to that of V; denote its 
configuration by A' . Next, melt the entire sample except 
for S, which remains frozen in configuration A. Now, 
allow the rest of the system to equilibrate; the region 
will now solidify into a new configuration B. If (with 
some probability) B ^ A' in some metrical sense, this 
means that the configuration in S determines (to some 
extent) that of a neighboring region proportional to 
its own size. Such a scenario is expected for quenches to 




FIG. 1: Sub-extensivity implies static correlation: The sur- 
roundings of identical patches are constrained. 

To motivate the definition of our correlation length, 
consider the following static variant of the preceding 
thought experiment. Beginning with an infinitely large 
configuration of a perfect equilibrium {i.e., ideal) glass, 
mark all regions of size £ isomorphic to S, that is, whose 
configuration is A; this is indicated shematically in Fig. 
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[T] (we consider only statistically isotropic systems here). 
We will loosely refer to all these regions as S's. The con- 
figurations of regions of volume aV contiguous with 
one of the E's is (with finite probability) determined by 
that of the neighboring S. Noting that S is represen- 
tative of regions of size it follows that if Q{V) is the 
number of possible configurations that a region of vol- 
ume V may take, then the constraint that S determines 
a finite fraction of XS' pairs may be expressed as: 

Q{aV)=-fn{V) . (1) 

If 7 is either constant or varies slowly with V, we have 
that 

n{v) - (2) 

with /i = -^j^- We note that for the case of a periodic 
arrangement, the state of is completely determined, 
so 7 = 1, implying ja = 0. On the other extreme, if 
the ideal glass state is random, then patches are inde- 
pendent and we have that \nQ(V) ex V: the complexity 
is proportional to the volume. Assuming that equal size 
patches occur with comparable probability, the distance 
V{i) we must go before we happen upon an identical 
patch of size V - £^ goes as V{£) - n^^^{V). In the 
example ^ above, V scales as power law V ^ i^^ to 
be contrasted with the uncorrelated case where it scales 
exponentially. We emphasize that the assumption that a 
region of volume V significantly restricts the possibilities 
of an adjacent region of size aV (with the concomitant 
power law for V) is rather strong. Although we regard 
this scenario as physically compelling, our definition of 
a diverging correlation length is more general, as it only 
requires that the number of embedded density profiles is 
subexponential in since then the repetition distance V 
will scale with £ subexponentially, which is all we need. 

With the previous discussion of pattern repetition in 
the ideal glass, we are now poised to define the coherence 
length. Given an (infinitely large) snapshot of a sam- 
ple, we pick a pattern of size £ and measure the spatial 
distance V{£) with which it appears to within a speci- 
fied precision, and repeat this procedure for the differ- 
ent patches of size £ (coincidence may be defined mod- 
ulo rotations, but this is inessential). Performing this 
procedure for different sizes, we plot, say, the average 
V{£) as a function of £. If the sample is above T^, or 
if it is below Tk but has not fully equilibrated to the 
ideal glass state, only small-scale patterns will repeat fre- 
quently, with large patterns recurring infrequently. That 
is, at some length £o there will be a crossover from subex- 
ponential behavior {e.g.; V ^ ^^), to a regime where the 
repetitions are much more sparse: V ~ e^^ . We define 
£o as the coherence length. 

We may relax the rigid requirement of strict identity 
and decide that two regions are congruent even if there is 
a number of mismatches. This allows us to work at finite 
temperatures, and to discuss correlation lengths in sys- 
tems undergoing general phase transitions. For example. 



in the case of a ferromagnet, the correct prescription is 
to identify two patterns if their overlap is larger or equal 
than the magnetization squared. One may also consider 
the generalization to coarse-grained, or even topological 
descriptions, which have recently proved interesting in 
the study of glassy systems [10 . 

It is important to note that sub-extensive complex- 
ity of embedded patches, such as that of Equation ([2|, 
does not imply that the infinite system is either periodic 
or quasiperiodic (these are characterized by a diffrac- 
tion spectrum consisting of (5— function Bragg peaks). 
To make this clear, we discuss a particularly well un- 
derstood class of systems, namely those which possess 
'recursive symmetry'. Examples of such systems are 
one-dimensional sequences of two symbols (1 and 0) 
which may be generated by means of a substitution rule, 
wherein a symbol is replaced by a specific word (com- 
posed of O's and I's) everywhere it occurs; generaliza- 
tion to higher dimensions is immediate. We consider 
three prototypes which characterize three types of spatial 
structure: 

Rule 1 : 1 ^ 10 ^ 10 (3) 

Starting from (say) 0, we generate sequences of length 2^ 
by iterating this rule n times; the result is clearly 

0101010101010101010101010101 

Continuing in this fashion, we generate a simple periodic 
sequence of arbitrary length. The second possibility is 
illustrated by 

Rule 2 : 1 ^ 10 0^1 (4) 

If we start from as before, we generate strings of length 
Fn^ where Fn is the n^^ Fibonacci number: 

1011010110110101101011011010110110 

This rule produces the Fibonacci sequence., which is 
quasiperiodic, whose Fourier transfom consists of a dense 
set of (5- functions. Finally, we consider 

Rule 3 : 1 ^ 10 ^ 01 (5) 

Iterating this rule n times produces a string of length 2^: 

01101001100101101001011001101001 

This sequence, known as the Thue- Morse sequence, has 
no ^-functions in its Fourier transform [TTl [121 US] • This 
kind of sequence also differs from quasiperiodic ones in 
that there are no finite translations that take a macro- 
scopic system into a similar configuration [14]. 

Consider the multiplicity of substrings of length £ of 
an infinite sequence generated by the above rules: in the 
first case there are only two possibilities, and in the last 
two 0{£) [T6l[T7j: all patches repeat in a distance V pro- 
portional to their size £. To see this, note that a subse- 
quence of size £ is completely included within the 'descen- 
dants' of a single site after k ex log(^) infiation steps; the 
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subsequence repeats in the descendant of a neighbouring 
ancestor - its 'cousin /c-times removed'. This argument 
is applicable to all systems obtained by substitutions. We 
conclude that the coherence length io is infinite for a per- 
fect sequence - even though there may be no Bragg peaks 
in the diffraction spectrum, so that by the conventional 
measure of scattering such systems would be considered 
amorphous. 

<©(£)> 



FIG. 2: The average distance V{i) between words of length i 
in a perfect Thue-Morse sequence of length 1000 (blue) and a 
sequence of length 1000 composed by concatenating strings of 
average length 50 randomly drawn from a perfect Thue-Morse 
sequence (red). 

Next, in analogy to a system which has not achieved 
the ideal glass state, we consider a sequence constructed 
by concatenating unrelated subsequences of varying 
lengths n, drawn from a distribution P(n), with aver- 
age n. It is easy to see that strings of length £ <^n will 
repeat often, because even in independent subsequences 
the same string will appear. In this case, V{i) will re- 
main close to its value in the ideal sequence. However, 
if ^ ^ n, the likelihood of overlap is exponentially small, 
by the central limit theorem. This is seen in Figure |2j 
where the average distance between strings of length i is 
plotted as a function of £ for a perfect and a randomly 
concatenated Thue-Morse sequence. 

These considerations may be generalized for higher di- 
mensions as well. Perhaps the best known 2D quasiperi- 
odic tiling is the Penrose tiling, but a system with similar 
properties is a Wang tiling [18] . In constructing a Wang 
tiling by successively laying tiles one next to another as 
a mason would, one is required to follow matching rules. 
One set of Wang tiles with their associated matching rules 
is shown in Figure [sja). An example of a system that is 
statistically isotropic is a tiling of identical triangles [20] 
called the 'Pinwheel' tiling. The Pinwheel tiling can be 
constructed by recursive substitution rules, so patterns of 
size £ repeat (modulo rotation) every P ~ £ in a perfect 
tiling, as for the ID sequences discussed above. Signifi- 
cantly, it is known that any packing which can be gen- 
erated by a recursive substitution rule admits a set of 
matching rules which only allows tilings consistent with 
the recursive symmetry [21; the set of matching rules for 
the Pinwheel tiling is known explicitly [20] . 

Having discussed the prefect order case £o = oo, we 
need to understand the nature of the defects that would 
reduce £0- Figure [4] illustrates how this comes about in 
the polycrystalline case; clearly £0 measures the typical 
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FIG. 3: (a) A set of Wang tiles with matching rules indicated 
by numbered edges. The rule for laying the tiles is that only 
edges of the same number may abut, (b) A state from a 
Monte- Carlo simulation of a Wang tiling quenched from high 
to low temperature: Only the 'hard' defects {i.e. mismatches 
of adjacent sides) remain; they are indicated by red lines. 




FIG. 4: Two patterns, one smaller and one larger than the 
crystallite length. The small pattern will recur with high fre- 
quency, but the large one will be exponentially rare, because 
it requires recurrence of the crystallite boundaries within a 
patch, which is very unlikely. 

crystallite size. While in this simple case the decoher- 
ence which decreases £0 is effected by well-defined grain 
boundaries, for general systems decoherence may occur 
due to less obvious defects. Consider for example the 
Wang tiling, which has the virtue that it can be expressed 
as a lattice model of spins |19. (taking 16 values), as- 
signing an energy corresponding to a mismatched edge - 
that is to say, the matching rules may be enforced by a 
suitable Hamiltonian. Consequently, such a system can 
be simulated at finite temperature (See Figure [sj^b)) to 
measure the growth of the coherence length. At finite 
temperature, the Wang tile system has energetic point 
defects - local violations of the matching rules. Start- 
ing from a high temperature state and annealing down 
to low temperatures, many easy-to-cure 'benign' defects 
disappear, leaving the system in a state with isolated 
'hard' defects (Figure |3jb)) that can only be removed by 
large-scale rearrangements of the order of the interdefect 
distance. Such behavior is also present in the Penrose 
tilings, where there are benign matching rule violations 
due to local tile flips, and non-trivial configurations of the 
phason field [22] which require large-scale rearrangements 
to eradicate. In this case, as we have verified numerically 
for the Wang tile system, the correlation length £0 is of 
the order of the distance between 'hard' defects. The 
conclusion we draw from this is that grain boundaries, 
playing both the role of borders of coherent regions and 
simultaneously being the place where energy is locally 
concentrated, is a peculiarity of crystals. Further work is 
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necessary to elucidate the various mechanisms for loss of 
coherence in quasicrystals, and nothing is known about 
such mechanisms for 'amorphous' recursive systems. 

It is worth relating the interesting concept of 'hyper- 
unifomity' (or 'superhomogeneity'), applied in the con- 
text of amorphous sphere packings [23l |24] , with the cur- 
rent work. A system is said to be hyperuniform if the 
square fluctuations in particle number between regions 
of identical volume V scales more slowly than V. One 
can restate this in terms of recurrence of patches of vol- 
ume V having the same particle number, while the cor- 
relation length defined here entails recurrence of actual 
patterns [24 . Note that a random stacking of hexagonal 
layers of hard spheres is hyperuniform, but at least in the 
direction normal to the layers, it has a finite io. 

The proposal that io diverges in the glass phase may 
seem to contradict the notion that glasses are 'spatially 
chaotic', since truly chaotic solutions do not have the 



frequent pattern repetitions we are claiming here. We 
argue that repetitions only appear when one considers 
the lowest free-energy states. A similar situation was 
already encountered in charge density waves [Ml [15], 
where chaos reduces to quasiperiodicity as soon as one 
restricts to the lowest energy solutions. 

To conclude, one may make the strong and intriguing 
conjecture that ideal glassy states have some recursive 
symmetry that can be constructed with some (yet un- 
known) infiation rule applying to continuous densities. 
In this case, it would be the emergence of this symmetry 
which is responsible for the diverging time scale. 
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